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Abstract. This paper is concerned with the stability of stationary solutions 
of the conservation law dtu + divyA{y, u) — AyU = 0, where the flux A is 
periodic with respect to its first variable. Essentially two kinds of asymptotic 
behaviours are studied here: the case when the equation is set on M, and 
the case when it is endowed with periodic boundary conditions. In the whole 
space case, we first prove the existence of viscous stationary shocks - also called 
standing shocks - which connect two different periodic stationary solutions to 
one another. We prove that standing shocks are stable in , provided the 
initial disturbance satisfies some appropriate boundedness conditions. We also 
extend this result to arbitrary initial data, but with some restrictions on the 
flux A. In the periodic case, we prove that periodic stationary solutions are 
always stable. The proof of this result relies on the derivation of uniform L°° 
bounds on the solution of the conservation law, and on sub- and super-solution 
techniques. 

Keywords. Viscous shocks; shock stability; viscous scalar conservation laws. 
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1. Introduction 

This paper is devoted to the analysis of the long-time behaviour of the solution 

u £ Ci[0,oo),Ll^{Q)) n L,^^([0,cx)),L°°(Q)) of the equation 

dtu + dYVyA{y, u) — AyU = 0, i > 0, y G Q, 
uit=o = "0 e L^{Q). 

Above, Q denotes either E or T^, the TV-dimensional torus (T^ = K^/[0,1)^), 
and A € Wl^^(j'^ X M)^ is an TV-dimensional flux (with = 1 when Q = M). 

Heuristically, it can be expected that the parabolicity of equation ^ may yield 
some compactness on the trajectory {u{t)}t>o- Hence, it is legitimate to conjecture 
that the family u{t) will converge as < ^ oo towards a stationary solution of |(T]). 
Such a result was proved when Q — by the author in [6] for a certain class 
of initial conditions, namely when uq is bounded from above and below by two 
stationary solutions of ^ . Such an assumption is in fact classical in the framework 
of conservation laws which admit a comparison principle: we refer for instance to 
[2j, where the authors study the long time behaviour of the fast diffusion equation, 
and assume that the initial data is bounded by two Barenblatt profiles. The same 
kind of assumption was made in the context of travelling waves by Stanley Osher 
and James Ralston in [17]; let us also mention the review paper by Denis Serre [19], 
which is devoted to the stability of shock profiles of scalar conservation laws, and 
in which the author assumes at first that the initial data is bounded from above 
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and below by shifted shock profiles. Nonetheless, in |[9j (see also [HITS]), Heinrich 
Freistiihler and Denis Serre remove this hypothesis, and prove that shock stability 
holds under a mere assumption on the initial data. 

The goal of this paper is to extend the result of [6] to arbitrary initial data, 
that is, to prove that solutions of ^ converge towards a stationary solution for 
any initial data uq S L°°(T^). We also tackle similar issues on the stability of 
stationary shock profiles in dimension one, when the equation is set on the whole 
space case {Q = R). Thus, a large part of the paper is devoted to the proof of 
the existence of shock profiles, and to the analysis of their properties. We will 
see that the question of shock stability reduces in fact to the stability of periodic 
stationary solutions of ^ in L^(R). Unfortunately, we have been unable to prove 
that periodic stationary solutions of Jl]) are stable in L^(R) for arbitrary fluxes. 
Hence we have left this issue open, and we hope to come back to it in a future 
paper. 

The proof of stability in the periodic setting relies strongly on the derivation of 
uniform L°° bounds on the family {u{t)}t>o. In the whole space case, the flrst step 
of the analysis is to prove the property for initial data which are bounded from 
above and below by viscous shocks; in fact, this result is similar to the one proved 
in |6j, and uses arguments from dynamical systems theory, following an idea by 
S. Osher and J. Ralston [17] (see also \W\ IT]). But the derivation of uniform L°° 
bounds is not sufficient to obtain a general stability result in the whole space case. 
Thus the idea is to use existing results on the long time behaviour of transport 
equations, which were obtained by Adrien Blanchet, Jean Dolbeault and Michal 
Kowalczyk in [5, and to apply those in the present context. 

Throughout the paper, we use the following notation: if u e L^(T^), 



Following [13], for a £ (0, 1), we define, if / is an interval in (0, oo) and is a 
domain in R^, 




We denote by 



(Q) the set of intergrable functions with zero mass: 




i7-'"(/ X n) 



{ueC{Ixn), ||m||^„/2.„(jxi^) < cx)}. 



where 



ff7'°(/xi2) 



max 

{t,x)£ixn 



\u{t,x)\ 




u{t',x)\ 



+ 



sup 
ix',t')eixn, 

\x—x'\<p 



|u(i, x) — u{t, x')\ 
\x - x'l" ' 



{x' ,t')eixQ, 
\t-t'\<p 



above, p is any positive number. We also set 




Eventually, for / e Lj"^ 



loc 



(R), heR,we set thJ = /(• + h). 
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2. Main results 

Before stating our main results, we recall general features of equation Jl]), to- 
gether with some facts related to the stationary solutions of this equation. 

In the rest of this paper, we denote by St the semi-group associated with equation 
([l]). Notice that St is always well-defined on L°°{Q), thanks to the papers by 
Kruzkov [Umi]- Moreover, we recall that the following properties hold true (these 
are called the Co-properties in |T9j): 

• Comparison: if a, 6 e L°°(Q) are such that a < b, then Sta < Stb for all 
t > 0. 

• Contraction: if a, 6 G are such that a— 6 G then 5fa— S'f6 e 
L^iQ) for alH > and 

\\Sta- Stb\\Li <\\a-b\\Li > 0. 

• Conservation: if a, 6 £ L°°{Q) are such that a — & e L^{Q), then Sta — 
Stb e L\Q) for alH > and 



(Sta - Stb) = (a^b) > 0. 
Jq 

Thanks to the Contraction property, the semi-group St can be extended on L°°{Q) + 
L^{Q). The so-called "Constant property" in |19] is not true in the present setting, 
since the flux A does not commute with translations. In other words, constants are 
not stationary solutions of equation H]) in general. The existence of periodic (in 
space) stationary solutions of (U) was proved by the author in [5], and we recall the 
corresponding result below: 



Proposition 2.1. Let A e W^^^iT^ x M)^. Assume that there exist Co > 0, 
m e [0, oo), n e [0, ^) when N>3, such that for all {y,p) e x M 

(2) \dMy,p)\<Ca{l + \pr) Vl<i<7V, 

(3) |div,^(y,p)| <Co(l + br). 
Assume as well that one of the following conditions holds: 
(4) 

/ fN + 2 \ 

m^OorO<n<lor ln< min ( — — — , 2 j and 3po £ K, divyA{-,po) = 



Then for all p E M., there exists a unique solution v{-,p) E Hpg^{T^) of the 



equation 

(5) - Ayv{y,p) + divyA{y,v{y,p)) = 0, {v{-,p))^p. 
The family (w(-,p))j,gR satisfies the following properties: 

(i) Regularity estimates: For all p E R, v{-,p) belongs to ^^^^^{T^) for all 1 < 
q < +00 and additionally 

(6) Vi?>0 3Cr>0 ypE[-R,R] \\vi;p)\\w2..(^T«j<CR. 
(ii) Growth property: if p > p' , then 

v(y,p) > v{y,p') VyET^. 
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(iii) Behaviour at infinity: assume that 
(7) snp\\dvA{-,v)\\j^^„M) < +00. 

Then 

lim sup v{y,p) = — cx), lim inf v{y,p) — +00. 

2.1. j4 priori bounds for solutions of scalar conservation laws. Our first 
result is concerned with the derivation of a priori bounds in L°° which are uniform 
in time. Notice that such a result is not trivial in general: in the homogeneous 
case, that is, when the flux A does not depend on the space variable x, this result 
follows from the comparison principle stated earlier. However, in the present case, 
this argument does not hold, since constants are not stationary solutions of H]). 
Of course, if there exists a constant C such that uq < v{-, C), then the comparison 
principle entails that StUg < v{-,C). Hence, the derivation of a priori bounds is 
easy when the initial data is bounded from above and below by solutions of equation 
^ . Consequently, the goal of this paragraph is to present similar results when the 
initial data does not satisfy such an assumption. 

Proposition 2.2. Assume that the flux A satisfies the assumptions of Proposition 
\2.1\ Assume also that for all K > 0, there exists a positive constant Ck, such that 
for all V e [-K, K], for all w eR, 

\divyA{y,v + w)-dWyA{y,v)\<CK{\w\ + \w\"), 

\d.Aiy,v + w)-d,Aiy,v)\<CKi\w\ + \wn, 

where n < {N + 2)/N. 

Let uo £ L°°{Q), and assume that there exists a stationary solution Uq G 
of (51) such that uoeUo+ L\Q). 

Then 

sup ||S'tUo||L~(Q) < +00. 
t>0 



Notice that in the above proposition, we do not assume that the stationary 
solution Uq is periodic. Thus Uq is not necessarily a solution of equation |[5]), and 
may be, for instance, a viscous shock proflle (see Proposition 12.41 below). In the 
periodic case, any function uq S is such that uq — w(-,0) S i^(T^), and thus 
the result holds for all functions in L°°. 

2.2. Stability of stationary periodic solutions in the periodic case. The 

derivation of uniform a priori bounds for the solutions of equation ^ allows us to 
extend the stability results proved in [6] to general classes of initial data. Let us 
flrst recall the stability result of [6]: 

Proposition 2.3. Assume that the flux A satisfies the assumptions of Proposition 
\2.1[ Let Uo G L°°(T^) such that there exists Pi, (32 & K such that 

(9) f(-,/3i)<uo<f(-,/32). 
Then as t 00 
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It was also proved in that under additional regularity assumptions on the flux 
A, the speed of convergence is exponential, due to a spectral gap result. 
We now remove assumption ([9]) thanks to Proposition 12.21 

Theorem 2.1. Assume that the flux A satisfies the assumptions of Proposition 
rO together with {S]). Then for all uq £ L°°(T^), ast-^oo, 

StUo^v{;{uo)) mL°°(T^). 

The proof of this result relies mainly on Proposition 12.21 and on sub- and super- 
solution methods based on the Comparison principle. Once again, it can be proved 
that the speed of convergence is exponential, provided the flux A is sufliciently 
smooth. For details regarding that point, we refer to [6]. 

2.3. Existence of viscous shocks. We now consider equation ^ set in Q = M. 
Our goal here is to prove the stability of a special class of stationary solutions, called 
"standing shocks". By analogy with the deflnition in [l9j of shocks in homogeneous 
conservation laws, a standing shock is a stationary solution U of equation ^ which 
is asymptotic to solutions of equation ([5]) at inflnity, namely 

3{pi,pr)eR^ lim iU{y)-viy,pi))=0, lim ({/(y) - = 0. 

y — * — oo y — ^+oo 

Because of the spatial dependence of the flux A, it does not seem to be possible 
to restrict the study of general shocks to standing shocks. For that matter, we wish 
to emphasize that the deflnition of a viscous shock with non-zero speed should not 
be exactly the same as in [l9|; indeed, it can be easily checked that if 

x) — U{x — st) 

is a solution of ([l]), then s = necessarily. Thus, for s ^ 0, a shock proflle is a 
solution of H]) of the form 

u{t, x) = U{t, X — st), 

where for all t, U{t) is asymptotic to solutions of equation ([5]) at inflnity. This 
is related (although not equivalent to) the deflnition of travehng pulsating fronts, 
see for instance the paper of Xue Xin [20] The existence of non-stationary shock 
proflles and their stability is beyond the scope of this paper, and thus, we will focus 
on standing shocks from now on, sometimes omitting the word "standing". 
Our flrst result is concerned with the existence of viscous shocks. 

Proposition 2.4 (Existence of stationary shock proflles). Assume that there exists 
p^ ,p^ G M such that p~ < p^ and 

(10) A{p+) = Aip-) := a, 

and define v± := v{-,p^). 
Let Uq G M such that 

v-{0) <ua< v+{Q), 
and Zei u : / — > R be the maximal solution of the differential equation 

(11) u (x) = A{x,u{x)) — a, 

(12) U|^=o = "0- 
Then u satisfies the following properties: 

(i) The function u is a global solution of JH]); in other words, / = R. 
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(ii) For all x&R, 

v-{x) < u{x) < v+{x); 
(iii) There exist qi, qr GM. such that qi ^ qr, A{qi) = A{qr) = a, and 
lim {u{x) — v{x, qi)) — 0, lim {u{x) ~ v{x,qr)) — 0. 

X — * — oo X — ^ + oo 

As a consequence, the solution u of (fTT |) - (fT2l) is a stationary viscous shock profile 
of equation ^ . 

Remark 1. (i) Assumption (fTO|) is the analogue of the Rankine-Hugoniot con- 
dition for homogeneous conservation laws. It is in fact a necessary condition, 
as demonstrated in Lemma [STTl below, 
(ii) In general, the asymptotic states v{qi),v{qr) are different from v{p'^),v{p^). 
Proposition 12.41 only ensures that 

A{qi)^A{qr)^A{p^). 

However, under an additional condition of Oleinik type, the asymptotic states 
can be identified: 

Corollary 2.1. Assume that the hypotheses of Proposition \2^\ are satisfied, and 
that the flux A is such that 

(13) Vpe A{p)^a. 
Then 

{quqr} = {p^,p^}. 

2.4. Stability of stationary shocks in the whole space case. We are now 

ready to state results on shock stability for equation Our first result is the 
analogue of Proposition l2.3l indeed, Theorem 12.21 below states that StUo converges 
towards a viscous shock, provided uq is bounded from above and below by the 
asymptotic states of the shock. In view of Theorem l2.1l it is natural to expect that 
this result remains true for arbitrary initial data. Unfortunately, we have not been 
able to prove this result in complete generality: we prove that stationary shocks are 
stable in provided stability holds (in -L^(IR)) for solutions of equation We 
also give explicit examples of fluxes for which the stability of shocks and periodic 
stationary solutions can be established. 

Theorem 2.2. Assume that the flux A satisfies the assumptions of Proposition 
\2.1[ Let pi,Pr S R such that pi ^ pr and A{pr) — A{pi) —: a, and assume that 
A,pi,pr satisfy Oleinik 's condition lfT3|) . 

Let U be a shock profile connecting v{pi) to v{pr). Let uq ^ U + L^(M) such that 
for almost every a; e M, 

(14) v{x,miii{pi,pr)) < uo{x) < u(x, max(p;,pr))- 

Then there exists a shock profile V connecting v{pi) to v{pr) and such that u S 
V + Ll{R). Moreover, 

^lim llS'tUo - "I^IIli(r) = 0. 

As outlined before, hypothesis lfT4l) should be compared with assumption 
Thus, the next step would be to prove that stability holds even when lfT4|) is false. 
In fact, we are only able to prove the following: 
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Proposition 2.5. Assume that the flux A satisfies the assumptions of Proposition 
\2.2l Letpi,Pr e K such that pi ^ Pr, A{pr) = A{pi), and such that (fT3|) is satisfied. 
Assume that the following assertion is true: 

(H) For p & {puPr}, for alluo &v{-,p) + Ll{M), lim ||S'tUo - i;(-,p)||ii(R) = 0. 

t — ^oo 

Let U he a shock profile connecting v{pi) to v{pr), and let uq £U + L^(M). Then 
there exists a shock profile V connecting v{pi) to v{pr) and such that u G V + L\(W). 
Moreover, 

^lim IjS'tUo - ^IIli(r) = 0. 

Remark 2. • As we will see in Section [3 hypothesis (H) can be relaxed into 

(H') For p e {pi,Pr}^ there exists 5 > Q such that for all uq G w(-,p) + io(K) , 
||uo - v{p)\\i <5^ lim llS'fUo - v{-,p)\\Lim) = 0. 

• In Section [71 we will prove the following result: for all p € M, there exists (5 > 
such that if uq € + io(^) satisfies ||mo — < (5, then 

lim ||5tUo - i'(-,p)||l~(r) = 0. 

Hence, in this case, for all e > there exists > such that St^u^ G v{-,p) + L\ 
and 

v{-,p-e) <uq < v{-,p + e). 
Consequently, hypothesis (H) can also be relaxed into 

(H") For all p eR, there exists S > such that for all uq G v{-,p) + Lq{R) 

v{-,p- S) <uo < v{-,p + d) ^ lim ||S'fUo - w(-,p)||li(R) 0. 

To sum up, denoting by (C) the conclusion of Proposition 12.51 (that is, shock 
stability), we have roughly 

(H) ^ (H") ^ (H') ^ (C). 

We now give an example when it is known that (H) is true. This example rehes 
on the analysis performed in the Hnear case by A. Blanchet, J. Dolbeault and M. 
Kowalczyk (see [3]). 

Proposition 2.6. Assume that the flux A satisfies the hypotheses of Proposition 
\2.2[ and let p E M. Assume that A is linear in a neighbourhood of v{-,p), i.e. 

36 e Ci(T^), > 0, Ve e {-V, v), Aiy, v{y,p) + ^ = A(y, v{y,p)) + b{y)C 

Then, provided a technical assumption on the moments of Stuo is satisfied (see 
(|42l) ), there exists S > such that for all uq G v{-,p) + Lq, 

Iko - v{-,p)\\i < ^ ^ lim \\StUQ - v{p)\\i - 0. 

t — *C30 

The assumption (|42| is a little technical to state at this stage, and is inherited 
from the analysis in [3j. However, as explained in [3j, this hypothesis is expected to 
be satisfied for a large class of initial data, so that in fact (|42|) is not a restriction. 

This allows us to give an expHcit case of fiux for which shock stability holds. 
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Corollary 2.2. Assume that the flux A is given by 

A{y,v)^V{y) + f{v), 
where V E C^(T^) and f E C^(R) is a convex function which is linear at infinity, 



I.e. 

\2 



f{v) ^ a+v ifv>A, 
f{v) — — a_u if V < —A. 



3(a_,a+) G [O.ooY, 3A>Q 

Then the following properties hold: 
(i) For all a > large enough, there exist {pi,pr) G such that pi ^ Pr and 
A{pi) = A{pr) = a, and there exists a shock profile U connecting v{pi) to 

V{pr). 

(ii) Let uq E U + L^(R). Then there exists a shock profile V such that u G V + Lq. 
Moreover, if l|42p holds for any vq G v{p^) + Lq{R), then limj^oo ||'5'tUo — 
^lli = 0. 

The plan of the paper is the following: given the similarity between the state- 
ments for periodic solutions when Q = T^, and stationary shocks when Q — M., 
we first prove the existence of standing shocks (i.e. Proposition 12. 4p and the shock 
stability result under boundedness conditions on the initial data (i.e. Theorem 12. 2p 
in sections [3] and |4] respectively. At this stage, we are able to treat both models 
simultaneously, and thus we prove Proposition l2.2l in SectionjH Section[6]is devoted 
to the proof of Theorem 12.11 and at last, we give in Section [7] further results on 
shock stability, including the proofs of Propositions 12.51 and 12.61 



3. Existence of one dimensional stationary viscous shocks 



This section is devoted to the proof of Proposition 12.41 together with a number 
of results related to viscous shocks which will be useful in the proof of Theorem 
12.21 These auxiliary results (monotonicity of shock profiles, integrability of the 
difference between two shock profiles, etc.) can be found in paragraph l3.3l We also 
give in paragraph l3.4l a few explicit examples in the case when the fiux A is convex. 

We begin with some comments on assumption IjlOp . 

3.1. Analysis of necessary conditions. 

Lemma 3.1. Let qi,qr G M, and let u G W^'°°{M) he such that 

u{x) — v{x, qi) as X ^ -co, 
u{x) — v{x, qr) ^ as X ^ +oo, 

^u" + ^{A{x,u{x)))^0. 
ax 

Then A{qr) = A{qi) =: a, and u satisfies 

—u + A{x, u{x)) = a. 

Proof. We deduce from the differential equation that there exists a constant C such 
that 

—u + A{x, u) = C, 

and the goal is to prove that A(q,.) = C = A{qi). We recall first that for all p G M, 
v{-,p) is a solution of 

-v'{x,p) + A{x,v{x,p)) = A{p). 
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Indeed, integrating ^ on R, we infer that for all p € M there exists a constant Cp 
such that 

Va; G M, ~v'{x,p) + A{x,v{x,p)) = Cp. 

Taking the average of the above equality over a period, we deduce that Cp — Mj>). 
As a consequence, we have 

(15) - (u(x) - v(x, qr)) + [A{x, u{x)) - A{x, v{x, qr))] = C - A{qr) 
ax 

Now, let (5 > arbitrary. There exists Xr > such that 

X > Xr =^ i\u{x) — v{x, qr)\ < S, \A{x, u{x)) — A{x, qr))\ < • 
Integrating (flSl) on the interval [xr , + 1] , we deduce that 

\C-A{qr)\ < 36. 

Since the above inequality is true for all 5 > 0, we infer that C = A{qr). The other 
equality is treated similarly. 

□ 

Remark 3. Notice that couples {pi,Pr) such that pi ^ Pr and A{pi) = A{pr) do 
not always exist. Indeed, consider the case of a linear flux A{x,v) — a{x)v, with 
a e C^(T). Then, for all p G M, we have v{x,p) = pm{x), where m is the unique 
probability measure on T satisfying 

-m/'{x) + -^{a{x)m{x)) ^0, xeT. 
ax 

The positivity of to is a consequence of the Krein-Rutman Theorem; we refer to [5] 
for more details. 

Therefore, for all p G M, 

A{p) = {av{-,p)) ^p{am) . 

Hence, as long as (am) ^ 0, A{p) ^ A{q) for all p, g G M such that p ^ q. In 
particular, if a is a non-zero constant, assumption (fTOl) is never satisfied. 

3.2. Proof of Proposition 12. 4L We begin with the a priori bound (ii), from 
which we deduce that u is a global solution. 

The inequality (ii) follows directly from classical results in differential equations; 
indeed, assume that there exists G / such that 

u{xo) > v+{xo); 

since u(0) < w+(0), there exists xi G [0,2:0] such that u{xi) = v+{xi). But u and 
v+ are solutions of the same differential equation, whence the Cauchy-Lipschitz 
Theorem implies that u = v+, which is false. Thus 

u{x) < v+{x) Vx G /. 

The lower bound is proved in the same way. 

As a consequence, we deduce that u remains bounded on its (maximal) interval of 
existence /. Using once again the Cauchy-Lipschitz Theorem, we infer that / = R, 
and thus m is a global solution. 

We now tackle the core of Proposition 12.41 First, since the flux A is T-periodic, 
the function u{- + 1) is also a solution of equation Hence the function x 1-^ 

u{x + 1) — u{x) keeps a constant sign on M, which entails in particular that for 
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all x G M, the sequences {u{x ± fc))fcgN are monotonous. Consider for instance the 
sequence of functions 

Uk : X G [0, 1] 1-^ u{x + k). 

According to the above remarks, the sequence [uk) is monotonous and bounded in 
L°°; hence for all x G [0, 1], Uk{x) has a finite Hmit, which we denote by Uoo{x). 
Moreover, thanks to the uniform bound (ii) and the differential equation ifTTj) . u 
belongs to 14^^'°° (K), and thus the sequence is uniformly bounded (with respect 
to fc) in M^^'°°([0, 1]). Thus u^o & W^'°°{[0,1]), and Uqo is a continuous function. 
According to Dini's Theorem, we eventually deduce that Uk converges towards Uqo 
in L°°([0, 1]). Notice that Uoc is periodic by definition, and passing to the limit 
in equation ((HJ, we deduce that Uoo is a solution of ifTT]) . Hence Uoo belongs to 
T4^i'°°(T) and satisfies 

-u'^ + ^{A{x,Uoo{x))) = 0, 
ax 

which means exactly that Uqo is a periodic solution of equation ([5]); according to 
Proposition 12.11 there exists p G R such that Uoo — v[-,p). Eventually, since Uoa is 
a solution of (fTTj) . we infer that a = A[p). To sum up, we have proved that there 
exists p G [p~,p'''], such that A{p) — A{p^)^ and such that 

lim sup \u{x + k) — v{x,p)\ = 0. 

The above convergence is strictly equivalent to u{x) — v{x,p) ^ as a; ^ oo, and 
thus the third point of the Proposition is proved. The limit as x ^ -co is treated 
similarly. 

3.3. Further results on viscous shocks. We have gathered in this paragraph 
some results which will be important in the proof of Theorem 12. 21 The first lemma 
gives a criterion allowing us to distinguish between the asymptotic states at ±oo. 

Lemma 3.2. Let pi,Pr G R such that A(j)i) — A{pj.), and let U he a shock profile 
such that 

lim \U{x) — v{x,pi)\ — lim \U{x) ~ v{x,Pr)] ~ Q- 

X — ^ — OO X — ^ + oo 

Then 

{d,A{;v{;Pl))) > 0, {d,A{;v{-,Pr))) < 0. 

Moreover, if one of the above inequalities is strict, then U converges exponentially 
fast toward the corresponding solution of equation 1^ ; for instance, if 

ar:= / dyA{y,v{y,pr)) dy < 0, 
Jt 

then for all a G (0, — a^), there exists a constant C'a such that for all y G [0, oo), 

\U{y) - v{y,pr)\ <Caexp{-ay). 

Proof. Since J7 is a shock profile and v{pi), v{pr) are solutions of equation ([5]), we 
have 

U'{x) ^ A{x,U{x)) 
dxv{x,pi) = A{x,v{x,pi)) - a, 

dxV{x,Pr) = A{x,v{x,Pr)) - Q!, 
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where a denotes the common value of A{pi) and A{pr). 

Consequently, the function U — v{pr), for instance, satisfies the linear equation 

(16) d^{U{x) - V(x,pr)) = b{x){U(x) - V{x,pr)), 

where 

h{x) = / dvA{x, tU{x) + (1 - T)v{x,pr)) dr. 
Jo 

Notice that since U converges towards v{pr) as a; — > +oo, we obtain 

(17) lim [h{x) ~ dyA{x,v{x,Pr))]=Q. 

On the other hand, equation (fT6|l impUes that 

U{x) - v{x,pr) = [C/(0) - v{0,pr)] exp (^^ b{y) dy^ . 
Once again, since U — v{pr) converges towards zero, we infer that 

(18) lim / b{y) dy ~ — oo. 

The first statement of the proposition follows easily from lfT7|) . (flSl) : indeed, assume 
that dr > 0. Then there exists a positive number K such that 

x> K ^ b{x) - dvA{x,v{x,pr)) > -y, 

and consequently, using the fact that x ^ dyA{x,v{x,pr)) is a periodic function, 
we obtain for x > K 

b{y)dy > / d^A{y,v{y,Pr)) dy - {x - K)^ 

K J K ^ 

> [x — K\ dr — x-^ — C 

> 4-a 

The above inequality is obviously in contradiction with (fTSj) . Hence < 0, which 
proves the first statement in the proposition. 

Now, assume that d^ < 0, and choose a e (0, — a^) arbitrary. As before, we pick 
K > such that 

X > K ^ b{x) — dvA{x, v{x,pr)) < —dr — a. 
We then obtain an inequality of the type 

b{y)dy < {-dr ~ a){x - K) + [x - K\dr + C 

K 

< -ax + C. 

Inserting this inequality back into the formula for U — v{pr) yields the exponential 
convergence result. 

□ 

The next result is concerned with the integrability of the difference between two 
shock profiles. 
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Lemma 3.3. Let pi,Pr € K such that pi ^ Pr and A{pi) — A{pr), and let U,V be 
two shock profiles with asymptotic states v{pi), v{pr). 
Then U-V e L\R). 

Proof. Set 

Uo := U{0), Vo F(0), 

and assume for instance that Uq < Vq. li Uq = Vq, then U — V according to 
the Cauchy-Lipschitz Theorem (see the proof of Proposition I2.4p . and the result is 
obvious. Thus we assume from now on that Uq <Vo. As a consequence, we have 

Vy e M, v{y,mm{pi,pr)) < U{y) < V{y) < v{y,ma.x{pi,pr)). 

We recall that the sequence {U{k))k£Z is monotonous, and converges towards 
v{0,pi) (resp. v{0,pr)) as k ^ —oo (resp. k +00). Hence, there exists an 
integer k^ £2, such that 

(19) C/o < < C/(fco), 

from which we infer that U <V <Tk„U . 

As a consequence, it is sufficient to prove that TkU— U is integrable, for all k £ Z. 

First, remember that TkU — U has a constant sign, since TkU and U are both 
shock profiles. Thus we only have to prove that the family 

{TkU - U) 



'-A 

remains bounded as A — > 00. A simple calculation leads to 

cA i-A i-A 

(TkU-U) = / U{y + k)dy- U{y)dy 

I J-A J-A 

U{y)dy- / U{y)dy 

k-A J-A 
k+A nk-A 

U{y) dy - / U{y) dy. 

I J -A 

Thus, recalling that i7 is a bounded function, we obtain 

i-A 



sup 

A>Q 



{TkU - U) 



<2k\\U\\ 



We deduce that TkU -U £ L^{M) for all A: e Z, and eventually that U-V £L^{m.) 
according to lfT9|) . 

□ 

The next result is in fact the first part of the statement of Theorem 



Lemma 3.4. Letpi,Pr £ K such that the assumptions of Theorem \2.S\ are satisfied, 
and let U be a viscous shock connecting v{pi) to v{pr). 

Let u £ U + L^ . Then there exists a unique shock profile V , with asymptotic 
states v{pi) and v{pr), and such that 

u£ V + Ll{W). 
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Proof. According to Lemma [3^ we already know that for every shock profile V, 
we have u — V^L^. Hence, the question is to find a shock profile V such that 



(20) / (u - y) = 0. 

Notice that such a shock profile is necessarily unique: indeed, the Cauchy-Lipschitz 
uniqueness principle entails that the difference of two shock profiles is a function 
which keeps a constant sign. Hence, if Fi, V2 are shock profiles satisfying j^iVi — 
V2) = 0, then Vi = V2. 

We now prove that there exists a shock profile V such that u — V ^ Lj(M). As 
before, we set p~ — min(p;,pr), = max(p/,pr)- For all ^ € (w(0,p^), w(0,p+)), 
we denote by the solution of 

V'{x)^A{x,V{x))-A{pi), 

Then, according to Proposition 12.41 and Lemma [3.3) for all ^, is a shock profile 
connecting v{pi) to v{pr), and additionally u — £ L^(R). Moreover, if ^ > 
then V^{x) > V^'{x) for all x; hence the function 



F : C G iv{0,p'),viO,p-^)) ^ / {u{x) - V^x)) dx 

Jr 

is well-defined and decreasing with respect to ^; using classical results on differential 
equations, it can easily be proved that F is continuous. We wish to find ^0 such 
that F{^o) = 0; thus it suffices to show that 

lim F(f) > and lim F(0 < 0. 

The above result is a direct consequence of Lebesgue's monotone convergence The- 
orem and of the fact that 

(21) VxeM, lim V^ix) ^v{x,p-). 

The same kind of result holds with w(p+). Indeed, let i? > be arbitrary, and let 
e > 0. Without loss of generality, assume that pr = p~ ■ Then there exists K & N 
such that 

X > K ^ v{x,Pr) < U{x) < v{x,Pr) + E. 

In particular, tk+ir\+iU is a shock profile which satisfies 

TK+yR\+iU{x) < v{x,pr) + £ Vx e [-R,R]. 

Let f := T/f_|_|^^j+i[/(0) = U{K + [R\ + 1). The Cauchy-Lipschitz Theorem entails 
that — T^-+iR\+iU . As a consequence, for all ^ < for all x € [—R, R], we have 

v{x,Pr) < V^{x) < V^{x) < v{x,Pr) + E. 

The convergence result l|2ip follows, and thus there exists a shock profile V such 
that UQ e y + Lj(R). 

□ 

The next lemma allows us to replace inequality (fT4l) by an inequality in which 
the upper and lower bounds are shock profiles, which will be useful in the proof of 
Theorem [221 in Section H 
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Lemma 3.5. Letpi,Pr such that the hypotheses of Theorem \2.2\ are satisfied. Let U 
he a shock profile connecting v{pi) to v{pr). Letu e L°°(R) such thatu e U + Lq(R) 
and assume that for almost every y G K, 

i;(y,min(pr,Pi)) < u{y) < v{y,ma.x{pi,pr)). 

Let £ > be arbitrary. Then there exists a function u'^ <E U + Lq{R), together 
with shock profiles Uj. connecting v{pi) to v{pr), such that 

Wu-u^Wli <£, VL-^u^ < C/+. 
Proof. First, since u — U e i^(K), there exists a positive number A'^ such that 

\u~U\ <£. 

x\>A' 

Hence, for > A'^ , we take = U{x). 

The definition of u"^ on the interval [—A^, A^] is sUghlty more technical, because 
of the various constraints bearing on u^. Once again, we assume that pi > pr in 
order to lighten the notation. We first consider a function £ C([— A'^, v4^]) which 
satisfies 



\u{x) — v^{x)\ dx < £ 



x\<A^ 



and such that 

v{x,Pr) < V''{x) < v{x,pi) Vx e [-A", A"]. 

We denote by a positive number such that 

v{x,pr) + < v^ix) < v{x,pi) - a" Vx e [-A^,A'']. 

Notice that can be chosen as small as desired. For further purposes, we choose 
so that 



{U~v{pr)). 



\<A^ 



The constraint u'^ £ U + Ll{M.) entails that the function should satisfy 



x\<A^ 

However, the function w"^ does not satisfy the above constraint in general: we merely 
have 



i 



{v^ - U) 



<A^ 



< 



< 



{v" - u) 



\<A^ 



(u-U) 



\<A' 



Iv" -u\ + 



x\<A^ 



x\>A^ 



< 2£. 



Assume for instance that /|^|<^e(w'^ — U) > 0. We then define a non-negative 
function p'^ G L°° {[- A"" , A^]) such that 



(22) 



t)^(x) - p^{x) > v{x,pr) + — a.e. on [-A", A" 



and 



/ {v- -p'-U) = 0. 

J\x\<A^ 
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Such a function p'^ exists provided 

a 



{V'~U)< Iv'- V{pr) 

\x\<A^ J\x\<A^ \ ^ 

and the above inequahty is equivalent to 



{U-V{pr)) > 
x\<A' ^ 

The previous condition is satisfied by definition of . Thus there exists a function 
which satisfies conditions (|22| . 
We then set 

u%x) = v%x) - p^x) for X e A""]. 
The construction is similar when /|a,|<^E(w'^ — U)<0. 

At this stage, we have defined a function u'^ + L}^ which satisfies 

vix,Pr) + Y <u'ix) <v{x,pi)-— Vxe [-A',A% 
u'{x)^U{x) Vx eR\[-A'',A'], 

and / |u — u^l < 4e. 



Now, by definition of the shock profile U, there exists a positive constant i?^ such 
that 

X> ^ \U{x) - v{x,pr)\ < —. 

Let fc+ be a positive integer such that /c+ > R^ + A^ . Then for all x e [— A"^, ^^], 
we have 

a" 

v{x,Pr) < Tk+U{x) < v{x,Pr) + ^ " '"^(^)- 

Similarly, there exists a negative integer such that for all x € [— ^^,74*^], 

a" 

u^x) < v{x,pi) ~ Y - '^i'-^^^^- 
Notice that tj,±(7 are also shock profiles. We now set 

Ul := sup(Tfc+C/, U), U'_ := M{Tk-U, U). 
Since shock profiles are ordered, the functions C/|. are viscous shocks, and 

C/1 < u= < Ul a.e. 

Hence the lemma is proved. □ 

3.4. An application: the convex case. This paragraph is devoted to the anal- 
ysis of specific examples for which the existence of shock profiles and their stability 
can be proved. 

Lemma 3.6. Assume that for all y G T, A{y, •) is a convex function. Then the 
homogenized flux A is convex. 

Furthermore, if A{y, •) is strictly convex for all y, then A is also strictly convex, 
and thus satisfies the Oleinik condition of Corollarv \2.1[ 

The convexity properties are proved in [15]. However, for the reader's conve- 
nience, we have reproduced the proof in Appendix B. Oleinik's condition is an 
immediate consequence of the strict convexity of the fiux A. 
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Example. Assume that the flux A is strictly convex in its second variable, and 
that the assumptions of Proposition 12 .41 are satisfled. Then, with the same notation 
as in Proposition 12.41 we have 

qi = p+ and qr =p~. 

Indeed, according to Corollary I2.H we have {qi,qr} — {p^,P~}- Moreover, since 
the flux A is strictly convex, dvA{y, •) is strictly increasing, and 

{d.Ai;vi;p-))) < {d.Ai;v{-,p+))) . 

Proposition 13.21 then allows us to conclude that p~ — qr , — qi- 

We now prove Corollary 12.21 (pending Proposition 12. 6p . Assume that the flux A 
is given by 

Aiy,p)=V{y) + fip), 

with V and / satisfying the assumptions of Corollarv l2.2l The existence of solutions 
of equation ^ follows immediately from Proposition 12. 1[ moreover, since the flux 
A is linear at infinity, hypothesis ([7]) is satisfied. As a consequence, for p > 
sufficiently large, we have 

A{y,v{y,p)) = V{y)+a+v{y,p) Vy G T^, 

and thus A{p) = {V) + a+p. Similarly, A{p) ~ {V) — a^p for p < with \p\ 
sufficiently large. These formulas entail that if a > is large enough, then, setting 
p± = ±(a — {V))/a±, we have A{p~) — A(p+) = a. Since A satisfies Oleinik's 
condition, we deduce that there exists a shock profile connecting v{p~) and f (p^). 

Additionally, if |p| is large enough, then / is linear, say, on the intervals [inf w(|p|)— 
l,oo) and (— cxd, supw(— |p|) + 1]. Thus, for all ^ € [—1,1], we have 

Myiv{y,p) + = V{y) + f{v{y,p) + = -4(?/,w(?/,p)) + sgn(p)asg„(p)C- 

Hence the fiux A satisfies the assumption of Proposition 12.61 for all p large enough. 
We infer that the solutions v{-,p'^) are stable by the semi-group St under small 
perturbations in Lq which satisfy (|42| . Point (ii) in Corollary 12.21 then follows 
from Proposition 12.51 and the remark following it. 

4. Stability of shock profiles in one space dimension - Part I 

This section is devoted to the proof of Theorem 12.21 Hence, throughout this 
section, we consider an initial data uq which satisfies lfT4|) . and such that uq G 
U + L^, where J7 is a stationary shock of equation H]). Using Lemma lS^ we deduce 
that there exists another shock profile V such that u G V + Then, using 

Lemma [331 together with the Contraction principle, we can restrict the analysis to 
the class of initial data uq such that 

(23) 3{U-,U+) shock profiles, U- < uq < U+. 

Indeed, assume that Theorem [521 holds for all vq e V+Lq such that (|23l) is satisfied. 
Consider now a function uq G V + Lq satisfying , and let e > be arbitrary. 
According to Lemma 13.51 there exists Uq G V + satisfying (|23l) and such that 
\\uo — UqWi < e. The contraction principle entails that for all i > 0, 

\\Stuo - V\\i < \\StUo ~ StulWi + WStu'o ~V\\i<e+ \\Stu'o - V\\i. 
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Notice also that by the Contraction principle, the function t ^ \\StUQ — V\\i is 
non-increasing, and thus has a finite limit as t — > oo. We infer that 

Ve > 0, lim \\StUo - V||i < e, 

t — >oo 

and thus SiUq converges toward y as t — > oo. 

There remains to prove Theorem 12.21 for initial data which satisfy l(23|) . As 
emphasized in Section [21 inequalities lfT4|) or l(23l) should be seen as the analogues 
of ([9]) in the context of shock stability. The proof of Theorem [221 in this case relies 
on arguments from dynamical systems theory, which are due to S. Osher and J. 
Ralston (see [l7j; similar ideas are developed by D. Amadori and D. Serre in p[]). 
The aim is to prove that the w-limit set of the trajectory Stuo is reduced to {V^}, 
by using a suitable Lyapunov function. Hence, we first prove that the w-limit set, 
denoted by il, is non-empty, then we state some properties of the w-limit set, and 
eventually we prove that fl = {V}. 

First step. Compactness in of the trajectories. 

Throughout this section, we set w{t) := Stuo — V. Notice first that by the 
comparison principle for equation ([1]), inequality i|23p is preserved by the semi- 
group St: for alH > 0, we have 

< Stua < U+. 

Hence, for alH > 0, 

U--U < wit) <U+-U. 

Since U+ ~ U and U- — U are integrable functions, the family {w{t)}t>o is equi- 
integrable in L^. Moreover, since U+ — U and U — C/_ are bounded, it follows that 
w is uniformly bounded in L°°. The function w satisfies a linear parabolic equation 
of the type 

dtw + dy{b{t,y)w) — dyyW — 0, t > 0, y £ R, 

with b e i°°([0,oo) X R). Theorem 10.1 in Chapter HI of [13] then implies that 
there exists a > such that for all tg > 1, for all R > 0, 

ll"(i)ll_ff"/2,»((to,to + l)x(-fl,_R)) < oo- 

Thus the family {w{t)}t>o is also equi-continuous in L^. 

Whence it follows from the Riesz-Frechet-Kolmogorov Theorem that the family 
{'w{t)}t>o is relatively compact in L^(R). Thus the w-limit set 

n:^\weV + L\R),3itn)neN:tn — > oo, St^uo Win L\R)} 

is non-empty. 

Second step. Properties of the uj-limit set fl. 

First, fl is forward and backward invariant by the semi-group St, meaning that 
for all t > 0, 

Stn = 

This important property is a generic one for cj-limit sets. It follows immediately, 
thanks to parabolic regularity, that all functions in fl are smooth: fl C H^^^(R), for 
instance. As a consequence, if W £ fl and wi{t) := StW, Theorem 6.1 in Chapter 
m of [I3] entails that wi G L^{[0,Tl H^{BR))nH^[0,T], L^{Br)) for all r,i?> 0. 

The second property which is important for our analysis is the LaSalle invariance 
principle (see [H]), which requires the existence of a Lyapunov function. In the case 
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of scalar conservation laws, a classical choice for a Lyapunov function is F[u] = 
\\u — V\\i. The Contraction principle entails that t F[StUQ] is non-increasing. 
Thus F takes a constant value on fl, which we denote by Cq. 

Eventually, using the conservation of mass, we deduce that is a subset of 

V + Ll 

Third step. Conclusion. 

We now prove, using the parabolic structure of equation (|24l) . that fl = {V}. 

Let Wo e be arbitrary, and let W{t) = StiWo). Notice that W{t) € for all 
t > 0, according to the previous step. Moreover, W satisfies 

dt{W ^V) + dy (A(y, W) - A(y, V)) - dyy{W ^V) = 0. 

Multiplying the above equation by sgn(T4^ ^ we obtain 

dt\W-V\+dy [sgn(W - V) {A{y, W{t)) - A{y, V))]-sgn{W -V)dyy{W ~V) = 0. 

Let (j) he a cut-off function, i.e. G C^(R), (j) > and = 1 in a neighbourhood 
of zero. For i? > 0, we set := (j){-/R). We now multiply the above equality by 
(t)R and integrate on [t,t'] x R. Recalling that J^\W{t) — V\ — Co for all t, we 
infer that for all t' > t> 0, there exists a function St.r ■ [0, oo) [0, oo) such that 

limij^oo £td'{R) = and 

rt' 

' sgn{W{s)-V)dyy{W{s,y)-V{y))^R{y)dsdy < et.t'{R). 



Thus, using a slightly modified version of Lemma [T] in the Appendix, we infer that 

Sgn(W(s) - V)dyy{W{s) -V)= dyy\{W{s) - V)\ 

almost everywhere and in the sense of distributions. Consequently, the function 
I W — y| is a non-negative solution of a parabolic equation of the type 

dt\W -V\+ dy{b{t, y)\W - V\) - dyy\W -V\=0, 

with b e L°°([0,ooM). We now conclude thanks to Harnack's inequality (see [8]): 
let G M be arbitrary, and let K be any compact set in R such that xa £ K. Then 
there exists a constant Ck such that 

|(Wo - V){xq)\ < sup \iWo - V)ix)\ < Ck inf |(W|,=i - V){x)\. 

xeK ^6-^ 

Now, (VF|s=i -V) <E Lin i/,^Q^(R), and thus there exists xi e R such that 

Wil,xi)-V{xi) = 0. 

Choose K such that xi € K. Then Wq — V vanishes uniformly on K, and in 
particular, {Wq — V){xq) = 0. Since xq was chosen arbitrarily, we deduce that 
Wo — V. Hence Q — {V}, and Theorem 12.21 is proved. 

5. Uniform in time a priori bounds for viscous scalar conservation 

LAWS 

This section is devoted to the proof of Proposition 12.21 As far as possible, we 
will treat both models simultaneously. We set 

w{t) :— StUQ — Uq, t >0. 

The function w satisfies the following equation 

(24) dtw{t, y) + d\YyB{y, w{t, y)) - l^yw{t, y)^0, t>0, yeQ, 
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where 

B{y, w) = A{y, Uo{y) + w) ~ A{y, Uo{y)), yeQ, weR. 

Due to the Contraction principle recalled in Section |2l it is known that w is 
bounded in L°°([0, oo), L^Q)), and 

(25) VteK+, \\w{t)\\Li <\\uq-Uo\\l^. 

The idea of this section is to use this uniform bound in order to derive uniform 
LP bounds on w for all p G [l,oo]. To that end, we proceed by induction on the 
exponent p. The first step is dedicated to the derivation of a differential inequality 
relating the derivative of the norm to a viscous dissipation term. The calculations 
are very similar to those developed in [5] to derive a priori bounds for solutions 
of equation Then, we use Poincare inequalities to control the norm by the 
dissipation. Eventually, we conclude thanks to a Gronwall type argument. 

Preliminary for the whole space case. 

We begin by recaUing some regularity results about the solutions of equation (U) 
in the case Q — R. According to the papers by Kruzkov [TTI, 112], it is kown that 
w e Li^^([0,oo),L°°(Q)). As a consequence, w e L^^{[0 , oo) , LP (Q)) for all p. 

Then, multiplying l(24|) by wx where x € Cq°(IR.) is an arbitrary non-negative 
cut-off function, and integrating in space and time, it is easily proved that for all 
T > 0, w satisfies an inequality of the type 



T 







\dyw{s,y)\'^xiy) dyds < Ct, 



where the constant Ct depends on T, HuiHiooQo.TjxR) and ||wt=o|ii; but not on x- 
We deduce that dyW e Lf^^{[0,oo), L'^{M.)). 

First step. A differential inequality. 

In this step, we treat the periodic and the full space models simultaneously; our 
goal is to prove an inequality of the type 



d 



V w ~ 



where 9 > 1 is arbitrary, n is the exponent appearing in ([8]), and the constants Cq 
and Cq depend on q, n, N, and ||C/o|| • 

To that end, we take q > 1, multiply l(24|) by and integrate on Q; we 

obtain 



1 



q + ldt Jq 



= g / Vyw{t,y) ■ B{y,w{t,y))\w{t,y)\''-Uy. 



Notice that all terms are well-defined thanks to the preliminary step. 

For (y, w) e Q X R, set 

bq{y,w)=q B{y,w')\wT^dw'; 
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then 



-q I Vyw{t,y) ■ B{y,w{t,y))\w{t,y)\''-Uy 
It 

divy {bq{y, wit, y))) + {d\Yybq){y, w{t, y))] dy 

w{t,y) 

{dWyB){y,w')\w'\'^~^ dw'. 

Thus, we now compute, for {y,w') G Q x R, 

diYyBiy,w') = div, [A{y, Uo{y) + w') - Aiy, Uo{y))] 

= (div^A)(y, Uo{y) + w') - {diYyA)iy, Uoiy)) 

+ VyUo ■ [{d,A){y, Uoiy) + w') - {d,A)iy, Uo{y))] . 

Consequently, according to hypothesis ([8]), we deduce that there exists a positive 
constant C depending only on ||?7o|!wi'°° and q such that 



Vyw{t, y) ■ B{y, wit, y))\w{t, y)\''-' dy 



< C 



9+1 



q+n 



Eventually, we infer that for all g > 1, there exist positive constants Cq, Cq such 
that for allt > 0, 



(26) I / \w{t)r^+cq I ^\w{t)\ — 



< C„ 



\w{t)\ 



9+1 



Second step. Control of norms by the dissipation term (Poincare inequalities). 

In this step, we treat the periodic case and the whole space case separately, and 
we begin with the periodic case. 

First, remember that for all p G (l,oo) such that j; > ^ — jj, there exists a 
positive constant Cp such that for all G H^^^{T^), 

(27) u-mp<Cp\\vn^. 



Taking ^ = [wl'a ^ we deduce that 

llu.||,<a 

where r G (l,oo) is such that 
(28) 



V k ~ 



+ ||w|| g + l 



1 1 

> 



r - q+l N{q + 1)' 

Now, the idea is to interpolate the and the L^+^ norms in the right-hand 

side of inequality (|26|) between and L*", where r satisfies the constraint above. 
It can be easily checked that when n < N + 2)/N, we have 

11 2 
n + q q+l N{q + 1)' 
hence the interpolation is always possible, and we have 

lkll«+i < ||w||}-"||HI", 

\\w\Ur.<Ml-^\\wf^, 
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where 



9- 



1 



a - 



a 
r ' 



1-/3 + 



Gathering all inequalities, we infer that 



II 119+1 
^ll-ll^+l 



Clklli' 



(g+n)(l-/3) 



V\w\ 



2 

2a 



CM': 



(9+l)(l-a)| 



,a(q+l) 



c\\wr; 



(9+n)(l-/3)| 



|/3(9+n) 



Remember that the norm is bounded. For the time being, we leave aside the 
norms of the right-hand side: those will be treated in the very last step. In 
order to control the right-hand side by the dissipation term in the left-hand side, it 
suffices to find r (and thus a and (3) such that 

2(3{q + n) 



(29) 



2a < 2, 



< 2. 



q + i 

Remembering the definition of we deduce that we have to find r e (q + l,oo) 
satisfying the two inequalities 

1 q+n-1 
1 > , 

r 9+1 

1^1 2 

r - q+1 N{q + 1)' 

This is possible if and only if the couple {n, q) satisfies 

g + n — 1 



1 



< 1, 
2 



< 1 



q + n- I 



q+l N{q + l) " q + l 

which amounts to the condition n < min(2, {N + 2)/N). In the case when = 1, 
this yields n < 2, which is more restrictive than the assumption of Proposition 12.21 
(n < 3). However, when = 1, the same arguments as in the whole space case can 
be used (see below), and lead to n < 3. Thus, under the hypotheses of Theorem 
12. H for a\\ q > 1, we may find r > max{q + 1, g + n) such that conditions (|28l) , ((29|l 
are fulfilled. The Cauchy-Schwarz inequality then implies that 
d I' - t „ii 2 



(30) 



dt 



\w{t)\ 



9+1 



V|u;(t)|- 



< 



c2(ii«;(t)rAi + ii^wrAi 

\ 2 2 



where the constant C2 depends on ||uo — ?7o||ii and the exponents pi,P2 on n, q 
and TV. According to the Poincare-Wirtinger inequality, we have 



V\w\ — 



> c 



I 1^ 
w 2 



|2±i 

w\ 2 



19+1 



IkllSl 
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Eventually, we deduce that for all g > 1, there exists constants Ci, C2,pi,P2 such 
that 

(31) I J |u;(t)|^+i + Ci 1 |u;(t)|^+i < {\\wit)r^ + \\wit)r^[ 

Let us now treat the one-dimensional model set in the whole space. In dimension 
one, the and norms control the norm. Hence we now interpolate the 

. . q+l 

two integrals in the right-hand side of l|26|) between and 



< \\w\\X\\wfS, 

2 

g+l g+271-1 



oo 



We use the following Poincare inequality, which involves the dissipation term in 
the right-hand side of (|26| (the proof of this inequality is classical and left to the 
reader: we refer to [TOj for the proof of similar inequalities): there exists a constant 
Cq, depending only on q, such that for all € n L°° n i?i(Mjl| 



< C„ 



1/3 



dy\w\ 



2\ 3(5 + 1) 



Consequently, there exist positive constants C,p such that for all w G n L°° n 



dy\w\ 



q+2n-l 
2\ 3(9+1) 



Hence, in order that the dissipation term controls the right-hand side of i(26|) . the 
exponent n should satisfy 

q + 2n-l 

which leads to the condition n < 3. Using Young's inequality, we conclude that 
(|30| is satisfied. Moreover, the Poincare inequality used above entails that for all 
A > 0, 

2\ 1/3 



< c\\w\\,4 



dy\w\ 



< 



dy\w\ 2 



Eventually, we deduce that inequality l(3T|) is also satisfied in the whole space case. 

Third step. Uniform bounds in for all q < oo. 

We now conclude thanks to Gronwall's lemma, using an inductive argument. 
Notice indeed that inequality l(3T|) implies that for all g > 1, 

(32) w G i°°([0, oo), L«(Q)) ^ w G i°°([0, oo), L29(Q)). 



^If w G L°° n H^{T), the corresponding inequality is 
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Indeed, assume that w G L°°([0, oo), L"?) for some q > I. According to l(3T1) . we 
have 

where the constant C2 depends on ||C/o|| and on ||ioo([o,oo),Li), so that, using 
Gronwall's lemma. 




Thus w e L°°([0,cx)),L29) and §^ is proved. Since w G L°°([0, 00), L^), we deduce 
that e i°°([0, 00), L?) for all q G [1, 00). 

Fourth step. Uniform bounds in L°° and W^'^. 

We now derive some L°° bounds thanks to parabolic regularity results. First, 
notice that in equation (|24l) . the flux B can be written as 

B{y,w{t,y)) = b{t,y)w, 

where ^ 

Kt^y)^ a{y,va{y) + Tw{t,y))dT. 
Jo 

According to the previous steps, b{t,y) e L°° {[0, 00), Ll^^{Q)) for all g > 0; in 
particular, in the whole space case, for all g > 1 there exists a constant Cq such 
that for all yo £ R, 

SUp\\b{t)\\L^y„^2,yo+2) < Cq. 
t>0 

We now use Theorem 8.1 in Chapter III of [l3j: we have, for all yo G Q, for all 

to > 1, 

where Qto,yo '■— (^0 ~ Ij^o + 1) (j/o ~ l^yo + 1) and q is some parameter chosen 
sufficiently large. The right-hand side is bounded uniformly in yo and to by a 
positive constant C, and we infer that for all yo & Q,to> 1, 

\w{to,yo)\ < C. 

Thus w e L°°([0,oo) X g). Using Theorem 10.1 in Chapter III of [13], we also 
deduce that there exists a > and a constant C > such that for all to > 1, for 
all xo e Q, 

lkllff°/2,= ((to,to + l)x(a:o-l,a;o + l)) - ^• 

As a consequence, we obtain 

lkllL~([l.oo)X"(Q)) < C. 

6. Long time behaviour of solutions for the periodic model 

Throughout this section, we assume that Q — , and we consider a solution 
u{t) = StUo of equation ^ {t > 0). Our goal is to prove, under the assumptions of 
Theorem l2.H that u{t) — v{{uo)) vanishes in L°° as t ^ 00. The idea is to prove in 
a first step the convergence for initial data which are bounded from above or from 
below by a solution of equation ^ , and then to extend this result to arbitrary initial 
data thanks to the L°° bounds proved in the previous section (see Proposition [221). 
We thus begin with the following Proposition: 
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Proposition 6.1. Let uq G L°°{T^) such that 

(33) 3poeR, uo{y)<v{y,pa) for a.e. y e T^" . 
Let u{t) — StUo for t > 0. Then, as t ^ oo, 

u{t)^v{-,{uo)) «nL°°(T^). 

Of course, the same result holds when the upper-bound is replaced by a lower- 
bound: 

Corollary 6.1. Let uq G L°°{T^) such that 

(34) 3^0 e M, uo{y) > v{y,po) for a.e. y G T^. 
Let u{t) — StUQ for t > 0. Then, as t ^ oo, 

u{t)-^v{-,{ua)) mL°°{T^). 

Proof of Proposition \6.1\ . According to the previous section (see Proposition 12. 2p , 

sup |ju(t)||ioo(T«) < +00. 

t>0 

Additionally, the Comparison principle yields 

u{t,y)<v{y,po) yt>0, Vy G T^. 

From now on, the proof is very close to that in [6], Section 2: we recall the main 
steps for the reader's convenience. Set 

U{t,y):^ sup uit',y), t>0,yGT^, 
t'>t 

p*{t) := M{p G R,v{y,p) > U{t,y) for a.e. y G T^} ,t > 0. 

Then U belongs to L°°([0, oo) x T^) (since u is uniformly bounded in time), and 
U is clearly a non-increasing function. Moreover, U satisfies 

U{t,y)<v{y,pa) Vt>0, Vy G T^. 

As a consequence, p*{t) is bounded from above by po, and p* is a non- increasing 
function. Moreover, p* is bounded from below, since for almost every y G T^, 

v{y,P*{t)) > U(t,y) > -||w||i^([o,oo)xT"), 

and thus 

yt > 0, p*{t) = {v{-,p*{t))) > -||w||ioo([o,oo)xT«)- 

Hence p* is a bounded decreasing function, and thus p*{t) has a finite limit, which 
we denote by p*, as i ^ oo. 

The idea is to prove that u{t) — v{-,p*) converges towards zero as t — > oo. Let 
e > be arbitrary. We first choose to > such that 

\\v{p*{t))~v{p*)\\oo<e yt>to, 

and then we pick p < p* and j/o G such that 

v{yo,P*) -£< v{yo,p) < U{to + l,yo) < v{yo,p*{to + 1)) < v{yo,p*) + e. 

Now, choose ti > to + 1 such that 

U{to + l,yo)-£ <u{ti,yo) < U{to + l,yo)- 
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By construction, the function 

V ■.is,y)e{-l,l)^T^ ^v{y,p*{to))-uih+s,y) 

is a non-negative solution of a linear diffusion equation of the type 

dtV + diyy{bV) ~ AyV = 

for some vector field b G L°°([— 1, 1] x T^)^. Hence by Harnack's inequality, there 
exists a constant C such that 

supv(-l-,y]<C mi V{0,y)<Ce. 

Thus, there exists a sequence of positive numbers (t„) such that lim„^oo tn = 
+00 and such that converges towards v{p*) in L°°. The contraction 

principle, together with parabolic regularity results, entails that the whole family 
u{t) converges. Eventually, we obtain that p* = (uq) by conservation of mass. 

□ 

The core of the proof of Theorem 12.11 then lies in the following argument: if 
uq e L°° is arbitrary, we set 

uo ■■= inf(uo,w(p)), 
u StUQ. 

The value of parameter p above is irrelevant. One can choose for instance p = 0, 
or p= (wo). 

The function uq obviously satisfies the assumptions of Proposition [Oj Hence as 

t oo, 

and thus there exists a positive time to such that for t >to, for all y G T^, 

u{t,y) > v{y, (uo) - 1). 

On the other hand, notice that uq < uq by definition, and thus by the comparison 
principle, 

u{t) < u{t) yt. 

Hence, for t> to, 

u{t) > V {{uo) - 1) . 

In particular, u{to) satisfies the assumptions of Corollarv l6.lt and thus, as t ^ oo, 

Stu{to) ^ V {{u{to))) . 

Since 

U{t) = St-toU{tQ) 

and {u{to)) — {uq) by the Conservation property, we deduce eventually that 

u{t) V {{uq)) as t ^ oo. 
Thus Theorem l2.1l is proved. 
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7. Stability of shock profiles in one space dimension - Part II 

This section is devoted to the proof of additional results on shock stability in the 
whole space case. We start by proving Proposition [231 a-nd then we prove that the 
conclusion of Proposition 12.51 still holds when (H) is replaced by (H'). We have 
not been able to prove that (H') is satisfied for arbitrary fluxes. Thus, at the end 
of this section, we prove Proposition 12.61 and thereby provide explicit examples for 
which (H') is satisfied. We also explain which difficulties are encountered when 
trying to prove (H'). 

We start by introducing some notation. Following [19], we denote by G the set 
of shock profiles connecting v{pi) to v{pr), and we set 

A:={ue (K), 3Ueg, ueU + L\R)} , 

Aq := {u e A, v{inm{pi,pr)) <u < w(max(p/,Pr))} • 

Our goal is to prove that for all uq G A, 

d{StUo,g) = 0, 

where d(u, A) denotes the distance from m to a set A. Notice that the Contraction 
principle easily entails that the function t i-^ d{StUQ,Q) is decreasing. Hence, its 
limit as < ^ oo exists; for all mq G A, set 

k){ua) := lim d{Stuo,g). 

t — ^oo 

Theorem 12.21 states that ioiu) — for all u £ Aq. Moreover, it follows from the 
Contraction principle that £o{uq) is a contraction, i.e. 

\io{u)-£oiv)\<\\u-v\\L^ Vu,weA 

Additionally, for alH > and for all u £ A, 

(oiu) = ioiStu). 

Similarly, we define, for all uq G A, 

ii{uo) := lim d{Stuo,Ao)- 

t — *oo 

The function £i is well-defined: indeed, the Comparison property entails that Ao 
is stable by the semi-group St- Consequently, by the Contraction principle, the 
function 1 1— > d{StUo,Ao) is decreasing and non-negative, and thus has a finite limit 
as t —>■ CO. Moreover, the functional £i enjoys the same properties as Iq: l\ is a 
contraction on A and l\{u) = l\{Stu) for all t > 0. Eventually, since Q C ^o, we 
deduce that 

4(u)<4(u) VueA 

7.1. Proof of Proposition [275l We now tackle the proof of Proposition l2.5l which 
is very similar to [19J, paragraph 3.5. Let wo € -4 be arbitrary. For all v G Aj^, we 
have 

4(wo) < + ||uo - ^'||l < ||uo - 

Thus for all uq £ A, 

loiuo) < d{uo,Ao). 

Replacing uq by StUQ in the previous inequality, we infer that for all uq G Aq, 

^(mo) < lim d{StUo,Aa) =^i(uo). 
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Thus io and £i take the same values on A, and it suffices to prove that 
(35) ei{uo) = lim d{StUo,Ao) = 0. 

t — >oo 

Notice that if u G A, then, with p+ = max(p/,pr), = niin(p;,pr), 

d{u,AQ)= {u-v{p+))^ ^+ {u-v{p-)) 

We now prove that assumption (H) implies ((35|) . According to Lemma lSTil there 
exists a shock proffie U such that u G U + Lq{R). We now define functions a+, a~ 
in v{p'^) + L\ and v{p^) + respectively, such that 

a~{y) < uo{y) < a+{y)- 

Let us explain for instance the construction of a+. If uo(y) > v{y,p'^), we set 

a+(2/) = uo(2/)- 

On the other hand, since u G U + and ?7 is asymptotic to v{p^),v{p~), we have 
- uo{y))+ dy> I {v{p+) -U)-\\uo- U\\i = +oo. 



Hence there is enough room, between the graphs of v{y,p^) and uo(2/) (restricted 
to the set where uo{y) < v(jj,p~^)), to insert a function 6+ such that 

uo{y) < v(y,p+) ^ uo{y) < b+{y) < v{y,p+), 

'i-uo<v(y,p+)iviy,P'^) - b'^iy)) dy^ I luo>"(y.p+)("o(2/) - v{y,p+)) dy. 

Jr 

On the set where uo{y) < v{y,p^), we define a+(j/) — h^{y). It is obvious that the 
function a+ belongs to w(p^) + L\ and that mq < a+. The function a~ is defined 
in a similar fashion. Thanks to the comparison principle, we have 

Sta^ < SfUo < Sta'^ Vi > 0. 

Consequently, 

(36) d{StUo,Ao) < \\Sta~^ - w(p+)||^i + \\Sta~ - v{p~)\\^^ . 

From the above inequality, it is clear that (H) entails l(35|) : if Sta^ — v{p^) vanish 
in L^, then ^i(uo) = 0. In other words, the stability of shock profiles follows from 
to the stability of solutions of equation ([5| in Lj. Thus, we now focus on the case 
when merely (H') is satisfied. 
Let S > 0. If Uo G ^ is such that 

||(uo - v{p+))+\\i < S, \\{uo - vip-))-\\i < 6, 

then by construction 

\\a+ - «(p+)||i < 25, - vip-)\\i < 26. 

And according to (H'), there exists So > such that if S < So, then 

lim \\Sta± - v{p^)\\i = 0, 

t — *oo 

and thus the right-hand side of ([36|) vanishes as i — > oo. Thus ^i(mo) = 0. 
Hence we now focus on the case where 

|l(uo - t'(p+))+l|i > So or \\{uo - vip-)).\\i > So. 
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We then define the function 

v{y,P^) + a+{uQ - v{p+)) if UQ{y) > v{y,p+), 

uo{y) ■■= { uo{y) if v{y,p-) < uo{y) < v{y,p+), 

v{y,P^) + a-{uo - v{p- j) ifuoiy) < v{y,p-), 



where 



^ r ii("o-"(f))^iu if _ „(p±))±||i > So, 

1 else. 



Since uq — uo G ^^(IR), uq G A. Moreover, 

\\uo-uq\\i = (1 - a+)||(wo - + (1 - a_)|l(uo - 

< d{uQ,Aa)-do. 
Notice that £i{uq) = 0. Since £i is a contraction, we have 

(37) ^i(mo) < £i{uo) + ||wo - wolli < d{uQ,Ao)-So. 
We now argue by contradiction. Assume that for all < > 0, 

\\{StUo - > So or \\{StUo - v{p~))^\\i > Sq. 

Then we may replace uo by StUo, for t > arbitrary, in inequality (|37|) . We obtain 

£i(mo) = £i{Stuo) < d{Stuo,Ao) - (Jo- 
Passing to the Hmit as t ^ oo, we infer 

^i(wo) < ^i(wo) - <5o, 
which is absurd. Hence there exists to > such that 

\\{StoUo ~ < So and \\{St„uo - v{p^))-\\i < Sq. 

We have already proved that £i{StoUo) = 0. We deduce that ^i(uo) = 0, and thus 
£{uo) = 0. 

Consequently, assumption (H') entails that £{u) = for all u ^ A. 

7.2. Stability of stationary periodic solutions in . We conclude this arti- 
cle by presenting some situations in which (H) or (H') hold true. We begin by 
explaining the linear case: assume that there exists a function b e Cper(IIi) such 
that 

Aiy,p)^biy)p V(y,p)e [0,1] xM. 

In this case, the stability of periodic solutions is a consequence of a result of Adrien 
Blanchet, Jean Dolbeault, and Michal Kowalczyk (see [Hill]): indeed, set a; = — (b), 
and let V G ^leii^) such that V' — {b) — b. Let p £ R be arbitrary, and let 
Uo e v{p) + Lq. Then, by Hnearity, w{t) ~: StUQ — v{p) solves an equation on the 
type 

dtw + dy{b{y)w) - dyyW = 0, 

t«lt=o = Wo e iJ(K). 

It is then easily checked that the function / defined by 

f{t, x) = w{t, X — Lot) 

solves 

(38) dtf{t, x) = d^^f{t, x) + 5, (^'(x - ujt)f{t, x)) . 
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This is precisely the case studied by Blanchet, Dolbeault and Kowalczyk. Let us 
recall briefly their method of analysis before stating their result. The flrst idea is to 
study the motion in the moving frame associated with the center of mass. Indeed, 
set 

x{t) := J xf{t,x)dx. 

Then it can be easily proved, using the linearity of the evolution equation and the 
periodicity of ^' , that 

lim -j-{t) = [ ^'m, 

t^OO dt Jq 

where m is the unique probability measure on [0, 1] solving 

(39) - m" + dx{{Lo + i)')m) = -to" + dx{bm) = 0. 

Set c :— {tp'm). The next idea is to perform a parabolic change of coordinates in 
the equation satisfied by /, in order to focus on the long-time behavior. Precisely, 
define U such that 

1 /, /- — X — ct 



Then U solves an equation of Fokker-Planck type, with a penalization growing ex- 
ponentially with time, and with coefficients which have fast oscillations for large 
times. Hence, this leads to the use of homogenization techniques, with the ad- 
ditional difficulty that the size of the oscillations in space depends on the time 
variable. An approximate solution is constructed thanks to a two-scale Ansatz. 
The convergence proof then relies on entropy dissipation methods. We are now 
ready to state their result. 

Proposition 7.1 (Blanchet, Dolbeault, Kowalczyk). Let wq G H L°°(K), and 
let f be the solution of l(38|) with initial data f\t=o — wq. Assume that there exists 
a constant Cq > such that 

(40) sup \ / \f{t,x)\{x~ct)Ux<Co. 

Then there exists a constant Ci, depending only on wq and Co, and a positive 
constant a, depending only on b, such that for all t >0, 

ll/Wlli<§- 



Remark 4. In fact, the result of Blanchet, Dolbeault and Kowalczyk is a little more 
accurate than the above proposition. Indeed, they prove that any non-negative 
solution / of (|38l) behaves asymptotically like 

/ /|t=o , , f x-ct \ 

, mix — Lot) /loo , , 

where ft-oo is a Gaussian function. In the present case, since wq & L^, the solutions 
and of (|38|) with initial data f^^Q = {wq)± have the same asymptotic 
behaviour, and thus w — f^ — f^ decays towards zero. 
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Consequently, in the linear case, assumption (H) is always satisfied, provided 
(|4Q1) holds. However, if the flux A is linear, standing shocks do not exist in general 
(see Remark [si . Thus, we now modify slightly the setting in order to use the 
results of the linear case, but in a non-linear context. Precisely, we now prove 
Proposition [2]6l which, as we have already stressed, provides an explicit example of 
shock stability without any assumption of the initial data except (|42l) (see Corollary 

E2D. 

Thus, let ^ be a non-linear flux which satisfles the assumptions of Proposition 
12.61 for some p S M. Let uq S v(j)) + Lq be arbitrary. Notice that we do not assume 
that 

\Wo - v{p)\\oo < V, 

so that the use of the linear setting is not straightforward. The idea is to prove, 
using dispersion inequalities, that 

(41) lim\\StUo-v{p)\\^^0, 

t — *C30 

provided ||uo — w(p)||i is sufficiently small. If the above convergence is true, there 
exists to > such that for t > to, for all y, {StUo){y) € [v{y,p) — r],v{y,p) +r]], and 
thus 

Aiy, Stuoiy)) = A{y, v{y,p)) + b{y)w{t, y) yt > h, Vy G M, 
with w{t) — StUQ — v{p). Consequently, for t > to, w solves a linear parabolic 
equation, and we can apply the previous analysis. The assumption on the moments 
of order four then becomes 

(42) 3^0 > 0, sup ^ [ \w{t, y)\ {y - jt)^ dy < oo, 

t>to + ^tj Jr 

where 7 := c — w. Notice that 7 = A'{p) in the present setting. Thus, we now focus 
on the proof of (|4T|) . 

As observed before, the function w{t) = StUQ — v{-,p) {t > 0) is a solution of 

dtw + dyB{y, w) - dyyW = 0, 

where the flux B is deflned by 

B{y,O^A{y,v{y,p)+0-A{y,v{y,p)), y e M, ^ e K- 

The idea is to linearize the flux B{-, ^) around ^ — 0, and to use energy methods. 
Let 

b{y) ^ d^A{y,v{y,p)), 
B{y,() = A{y,v{y,p)+0 - A{y,v{y,p)) - h{y)£,. 

Since A € W'^'°°(T x E), the flux B is quadratic in a neighbourhood of ^ = 0. 
According to Proposition 12.21 w is bounded in L°°([0, 00) x R), and thus there 
exists a constant C such that 

\Biy,wit,y))\<C\wit,y)\^ > 0, Vy G M. 

The function w solves 

(43) dtW + dy{b{y)w) - dyyW = -dyB{y,w). 

Following an idea of Philippe Michel, Stephane Mischler and Benoit Perthame (see 
[T6]l. we consider the invariant measure to, deflned by l(39|) . Notice that m = dv/dp 
in the present case, and there exists a positive constant C > 1 such that 

< TO < C, \dym\ < C a.e. 
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Moreover, according to jT6], the following identity holds 



dt m 



dy \ m 



= -2m 



(-) 



2 W ~ 

-2-dyB{y,w) 



Integrating the above equation on R, we obtain 



1 d 
2dt 



m 



+ m 



< C \w\ 



w 



< c 


w 


2 












m 


L4(m(y)dy) 


m 



L^{m{y)dy) 

Notice that for all p G [l,oo), the norm is equivalent to the U'(m{y)dy) norm. 
We now use the following Poincare inequality: there exists a positive constant C, 
such that for all (}> G ^^(R) n -ff^(R), there holds 



< c'II'/''IIl2^(b 



1/2 



(44) ||<^||l4, 

Taking (j> — w/m, we are led to 

1 d f w 2 r „ w 2 „„ I 

2 dt m m 

< c||woIIli( 

Now, if ||wo||li(b) is sufficiently small, we obtain 



m 
m 



(45) 



d_ 

dt 





w 


2 






/ m 

m 


m 




m 


m 



< 0. 



We then proceed as in pjj (Paragraph 1.1): using the Nash inequality together 
with the decay of the norm, we deduce that for all t > 



w{t) 



<C\\w{t)\\%' 



and thus we infer 



,wit) 



1/3 



<C\\wo\\l{' 



L2(m) 



,w{t) 



1/3 



d 

di 





w 


2 






w 




/ m 

m 


m 


+ 


ll^oiii. ( 


/ m 

Jr 


m 






< 0. 



Integrating the above differential inequality, we obtain eventually 

||u'(i)||L2(R) < C 



Thus the norm decays with an algebraic rate. 

We now use a parabolic regularity result, from which the decay of the L°° norm 
immediately follows. The key point Hes in the following inequality: there exists a 
constant C such that for alH > 1, 

\\w{t)\\^ < C\\w{t - l)h. 

Indeed, w satisfies an equation of the type 

(46) dtw + dy{a{t, y)w{t, y)) - dyyW = 0, 
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where the coefficient a is bounded in L°°{[0,oo] x M). Using an energy estimate, 
it can be easily proved that there exists a positive constant a, depending only on 
||a||oo, such that if W is any solution of (|46l) . then for any t > s > 

\\W{t)\\2<e^('~^^\\Wis)h. 

Moreover, according to Harnack's inequality, there exists a constant C such that 
for any non-negative solution W of l(46|) . for alH > 0, 



W{t,y)<C inf W{t + l,z)dz. 

ze[y-l.y+l] 



And if t > 1, 



/I \^ 
inf W(t+l,z)dz < -/ W'^(t + l,z)dz] 

z£ly-l,y+l] \^ Jy-1 J 

< Ci W^{t~l,z)dz\ 

< C\\Wit-l)\\L2(M)- 

Now, let t > 1 be arbitrary, and let := {w{t - 1))+, := {wit - 1))_. For 
s >t—l, consider the solution Wi of l|46p such that Wi\s=t-i — . The functions 
Wi are non-negative by the maximum principle. Consequently, for all y € R, we 
have 

W.{t,y)<C\\W^h. 
Since w = W\ — W2, we deduce that 

sup\wit,y)\<C\\w{t-l)\\2. 

The decay of the L°° norm follows. This concludes the proof of Proposition [2]6l □ 

In fact, the decay of all norms for p e (l,oo] is a general property, which 
is true even when the flux A does not satisfy the assumptions of Proposition 12.61 
However, if the flux A is not linear in a neighbourhood of v{p^), then we are unable 
to conclude to the stability of periodic solutions. Let us now explain briefly where 
the difficulty lies: a natural idea would be to treat the term i? as a perturbation in 
(|43l) , and to write a Duhamel formula of the type 



w{t)^S^wo+ / St, dyBi-,w{s)) 







ds. 



where is the (linear) semi-group associated with the equation 

dtW + dy{hw) — dyyW = 0. 

This is exactly the method used in [19] in order to prove the stability of constants 
in the viscous model. Thanks to the results of [3], it is already known that SfWQ 
decays in as i ^ 00. However, the method used in [T9j cannot be used here, 
essentially because of the complicated dependance of the constant Ci appearing in 
Proposition 17. II on the function wq. Indeed, looking carefully at the proof in [3], it 
can be checked that 

X 1/2 



^ "^0 111 \\\w^\\ihooJ Jr\\wo\\i V||wo||i^o ' 
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where hoc is a normalized gaussian function, Wq,Wq are the positive and negative 
real parts of wq, and 

Co ■■= sup ^ / \StWo\ix-A'ip)t)\ 

The dependance of Ci with respect to Co and the relative entropies of Wq and Wq 
is disastrous for the use of the Duhamel formula: indeed, one has to control, for 
instance, 

r idyBi;w{s)))+ ( dyB{;wis)))+ \ 

Jr \\dyB{-,w{s)))+\\l '^\\\dyB{-,w{s)))+\\lhoc) ' 

Consequently, another approach must certainly be chosen in order to prove the 
stability of periodic solutions in the general case. Given the complexity of the 
proof in the mere linear setting (see [3]), this question goes beyong the scope of 
this article. Also, we emphasize that it is not clear that the methods of [3] can be 
adapted to a nonlinear setting: indeed, the proof of convergence relies on the use of 
entropy dissipation techniques, which are more adapted to the Hnear case. We refer 
to [7] for additional results and techniques concerning the asymptotic behaviour of 
non linear viscous conservation laws in the homogeneous case. 

Appendix A - Proof of Lemma EH] 

Assume that the flux A is convex, and let pi,P2 G such that pi 7^ P2 , and let 
A G (0, 1). In the following, we set 

Vt{y) =v{y,pi), i = l,2, 
w — Xvi + (1 - X)V2, p = Api + (1 — A)p2, 
My) = v{y, Xpi + (1 - A)p2)- 
By definition of v{-,p) and of the homogenized fiux A, we have 

-vl + A{y,v,{y))^A{p,), 
-u' + A{y, u{y)) = A{\pi + (1 - A)p2). 
Consequently, using the convexity of the fiux A, we deduce that for all y e T^, 
(47) -w'{y) + A{y,w{y)) < -w'iy) + \Aiy,vi{y)) + {1 - \)A{y,V2iy)) 

- AA(pi) + (l- A)A(p2). 
Assume that A{Xpi + (1 — X)p2) > AA(pi) + (1 — X)A{p2), and write u, w as 

u = p + ,f\ w^p + g', 

with f,g£ Cpgj.(T^). Since / and g are defined up to the addition of constants, we 
can assume that f < g almost everywhere. Moreover, notice that 

sup {-g"{y) + Aiy,p + g' {y))) < inf (-/"(y) + A(y,p + /'(y))) . 

Thus there exists a > such that 

-5" + Aiy,p + g'{y)) + ag < -f" + A{y,p + f {y)) + af. 
Hence, by the maximum principle, we infer that g < f, which is absurd. Thus 

A{Xpi + (1 - X)p2) < XA{pi) + (1 - X)A{p2). 
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If the flux A is strictly convex, then inequahty l(47|) is strict for all y G 
(remember that the family v{y,p) is strictly increasing with p for all y G T^). 
Consequently, the same argument as above leads to 

A{Xpi + (1 - X)p2) < AA(pi) + (1 - X)A{p2) 
Appendix B 

1. Let w € L'^n L°°(R) such that w' e L'^{R) and w" € i^cW- 

Assume 

that w is such that 



lim / sgn{w{y))w"{y)(f>(^) dy = 
for all (j) € C^(R) such that cf) = I in a neighbourhood of zero. Then 



lim ^ I k'pl|u,|<5 = 0. 



1 

As a consequence, 

dyy\w\ = sgn{w)w" in V 



Proof. For 5 > 0, let 



sgn(x) if |a;| > S, 

X 

1 



^ else. 



Then ^ 
and for all i? > 0, we have, using the chain rule 



where <j)i{. = <j)[-/R). 

Since w' G LP' , we infer 



Thus the above term vanishes as i? oo, uniformly in S. 
On the other hand. 



L2. 



lini J w"tps{w)(j)ii = J w"sgn{w)(f>R, 
and the right-hand side vanishes as i? ^ cx) by assumption. We deduce that 

lim limsup / \w'\'^ip'g{w)(j)ii = 0. 

Now, since the integral / \w'\'^'ip'g{w)(j)ii is non-negative and increasing with re- 
spect to R, we deduce that 

lim / \w'\^^'s{w)cj)R = Vi?, 



and thus the first part of the lemma is proved. 

^^^^ 



Consider Ss € W^^^{M.) such that 



S's = and ^^(O) = 0, 
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where the function i/'cS was defined earlier. Then 



and according to the chain rule, 



dyySsiw) 



w"i1js{w) + \u)'\ 




6 



Passing to the limit in the sense of distributions in the above equality yields 



I am very grateful to Denis Serre, for encouraging me to work on these questions 
in the first place, and for very interesting and fruitful discussions. I also wish to 
thank Jean Dolbeault, for his helpful insight of the long-time behaviour of linear 
diffusion equations. 
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